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Abstract 

Experimental and modeling/simulation studies of phase equilibrium and growth morphologies of novel 
polymer-dispersed liquid crystal (PDLC) mixtures of PS (polystyrene) and liquid crystals that ex- 
hibit a direct isotropic/smectic-A (lamellar) mesophase transition were performed for PS/10CB (decyl- 
cyanobiphenyl) and PS/12CB (dodecyl-cyanobiphenyl) . Partial phase diagrams were determined using 
polarized optical microscopy (POM) and differential scanning calorimetry (DSC) for different composi- 
tions of both materials, determining both phase separation (liquid/liquid demixing) and phase order- 
ing (isotropic/smectic-A transition) temperatures. The Flory-Huggins theory of isotropic mixing and 
Maier-Saupe-McMillan theory for smectic-A liquid crystalline ordering were used to computationally de- 
termine phase diagrams for both systems, showing good agreement with the experimental results. In 
addition to thermodynamic observations, growth morphology relations were found depending on phase 
transition sequence, quench rate, and material composition. Three stages of liquid crystal-rich domain 
growth morphology were observed: spherical macroscale domain growth ("stage I"), highly anisotropic 
domain growth ("stage IP), and sub-micron spheroid domain growth ("stage III"). Nano-scale structure 
of spheroidal and spherocylindrical morphologies were then determined via two-dimensional simulation of 
a high-order Landau-de Gennes model. Morphologies observed during stage II growth are typical of di- 
rect isotropic/smectic-A phase transitions, such as highly anisotropic "batonnets" and filaments. These 
morphologies, which are found to be persistent in direct isotropic/smectic-A PDLCs, could provide new 
functionality and applications for these functional materials. 



1 Introduction 

Polymer dispersed liquid crystals (PDLCs) are func- 
tional composite materials consisting of an isotropic 

* Institute of Materials Science and Technology (INTEMA), 
University of Mar del Plata and National Research Council 
(CONICET), J.B. Justo 4302, 7600 Mar del Plata, Argentina 

t Department of Chemical Engineering, McGill University, 
Montreal, Quebec H3A 2B2, Canada 



polymer matrix and a liquid crystalline droplet phase 
[1-3]. These functional materials are used in states 
where liquid crystalline ordering drives the primary 
function of the material as in switchable windows, 
displays, spatial light modulators, tunable filters and 
other devices [1]. PDLC materials can also be ther- 
mally actuated to transition from opaque to trans- 
parent states, allowing for other types of applica- 
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tions such as thermal sensors and active elements in 
thcrmo-optical memory devices [4]. 

The formation of PDLC materials can be driven by 
three phase separation mechanisms: polymerization- 
induced (PIPS), thermally-induced (TIPS), and 
solvent-induced (SIPS). Phase transition phenom- 
ena exhibited by PDLCs are complex, since a com- 
bination of phase separation (chemical de-mixing), 
polymerization (as in PIPS), and phase ordering 
(liquid crystallinity) can be involved. Thus PDLC 
phase transformation processes involve coupled con- 
served (concentration) and non-conserved (liquid 
crystalline) order parameters. Phase transitions in- 
volving multiple order parameters have been shown 
to lead to unique material architectures, such as col- 
loidal crystals dispersed in liquid crystal matrices [5] . 

Under isothermal conditions and in absence of con- 
vective flow and polymerization (as in TIPS), the ki- 
netics of the coupled order parameter process involv- 
ing demixing and phase ordering is described by: 
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where C is the polymer concentration, Mj are mobil- 
ities, F is the free energy, and is the order param- 
eter of the mesophase. Since these order parameters 
(C, 0) follow different kinetic laws, their interaction 
results in a rich and diverse set of anisotropic growth 
morphologies (below the liquid crystalline transition 
temperature), defect formation, and domain textur- 
ing (gradients in liquid crystal order). As demon- 
strated in ref. 6, intrinsic coupling terms between or- 
der parameters in the free energy F of the form: 



Fee = icV6 • WC 



(3) 



contribute to the shaping of anisotropic domains, as 
well as texturing (Lc is a coupling parameter). In ad- 
dition, the transformation growth laws, usually given 
in terms of power laws R cx t n (typically n = ^ for 
conserved and n = 1 for non-conserved under deep 
undercooling [7,8]), will exhibit non-trivial dynamics 
due to the coupled order parameters. 

The majority of research in the field of PDLCs has 
focused on liquid crystals which exhibit partial ori- 
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Figure 1: Schematics of the (a) isotropic (disordered), 
(b) nematic, and (c) smectic-A phases. 

entational order [4,6,9-15], the nematics (see[TJi-b). 
Higher-order liquid crystals also exist, such as smectic 
and columnar mesophases, where both orientational 
and translational is present (see [Ik) . The thermo- 
physical characterization of these higher-order liquid 
crystals is less complete compared to nematic liq- 
uid crystals. This situation is amplified when con- 
sidering PDLCs involving smectic mesogens [16-18] 
in that eqns. [T][2]have to include positional and ori- 
entational order. Additionally, when considering the 
formation of smectic-A mesophases, it is found that 
some compounds exhibit a direct isotropic/smectic- 
A (JTJa, — >[lp) transition upon quenching, while others 
exhibit a sequential symmetry breaking cascade ([lji 
— > IIJ) — > uL). The direct isotropic/smectic-A tran- 
sition exhibited by both 10CB and 12CB and is a 
simultaneous symmetry breaking transition. It is of 
particular interest due to the complex growth mor- 
phologies [19, 20] (rods, filaments, "batonnet" struc- 
tures) and dual order parameter transformation ki- 
netics [21] that have been observed in past studies. 
When considering TIPS-based PDLCs, it is expected 
that the additional presence of mass diffusion will ex- 
pand the phenomenology already found when using 
nematogens [22]. In particular, there will be three 
binary coupling free energies (as eqn. [3] one con- 
served and two non-conserved), assuming isothermal 
conditions. 

PIPS has been the standard approach to produc- 
ing PDLCs in that the final PDLC material has a 
structure which is more thermodynamically stable 
when exposed to temperature changes and solvents 
[4,6,9-15,23]. In this work we use SIPS to gener- 
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ate the initial material, and reset the microstructure 
using TIPS. 

The overall objective of this work is to character- 
ize the thermodynamic and material transformation 
processes of novel polymer-dispersed liquid crystal 
(PDLC) mixtures of PS (polystyrene) and liquid crys- 
tals that exhibit a direct transition from a disordered 
to the lamellar smectic-A state. A comprehensive 
approach is taken that includes both experimental 
and computational techniques in order to access the 
multiple scales involved in PDLC phenomena. The 
specific objectives of this work on smectic PDLCs 
are divided into two areas: (1) equilibrium thermo- 
dynamics and (2) material transformation processes. 

The first objective, studying equilibrium thermo- 
dynamic properties of these novel smectic PDLC ma- 
terials, is accomplished by first experimentally de- 
termining a partial phase diagram. This includes 
determination of the isotropic/isotropic phase sep- 
aration temperature (Tjj) and isotropic/smectic-A 
phase transition temperature (Tat), or at high con- 
centrations of polymer the temperature at which both 
phase separation and ordering take place simultane- 
ously (Tj/ai). Following this, experimental results 
are then compared to computationally determined 
phase diagrams using the Flory-Huggins theory of 
isotropic mixing [24] and Maier-Saupe-McMillan the- 
ory [25], for both material systems. 

The second objective involves determining phase 
transformation growth morphology trends based 
upon thermodynamic pathway, guided by the phase 
diagrams determined in Objective 1. A general the- 
ory identifying three distinct growth regimes (stages 
I, II, and III) is used based upon research in PIPS 
nematic PDLC systems [26]. Two-dimensional sim- 
ulation and scaling theories are then used, based on 
experimental morphology observations, to shed light 
on the general nano-scale structure of spheroidal and 
sphero-cylindrical domains. 

The paper is organized as follows: the experimental 
procedure, models, and simulation conditions are pre- 
sented, followed by the quantitative thermodynamic 
results (both experimental and computational) . This 
is followed by a brief discussion of the experimen- 
tally observed phase transition behavior, referring to 
related past work studying lyotropic liquid crystals. 



Qualitative trends of morphological growth are then 
given with supporting POM images. Following dis- 
cussion of the macro-scale experimental observations, 
the nano-scale simulation results are presented. Fi- 
nally, conclusions are made, including guidance for 
future work. 

2 Experimental 

2.1 Materials 

Mono-disperse polystyrene (PS, Polymer Source), 
with molar mass M n = 20800 and M w = 22200 (pro- 
vided by manufacturer), was used as received. The 
glass transition temperature is 99°C, and the mass 
density at room temperature is 1.047 g/cm 3 . The 
LCs 4-decyl-4-cyanobiphenyl (10CB, Synthon) and 4- 
dodecyl-4-cyanobiphenyl (12CB, Synthon), were used 
as received. Both LCs show a isotropic / smectic tran- 
sition (no stable nematic phase) at 52.4 °C (10CB) 
and 59.6 °C (12CB). The melting temperatures are 
46.5 °C (10CB) and 48.2°C (12CB). 

2.2 Sample preparation 

The mixtures, with different compositions of PS and 
liquid crystal, were dissolved in THF and stirred 
manually until complete dissolution. For microscopy 
measurements, a small amount of the mixture was 
cast on a clean glass slide and the samples were left 
for 24 hours to allow for complete evaporation of the 
solvent (verified by tracking the weight of the sam- 
ple). Another glass slide was put on top of the first 
one, and the dry sample was sandwiched between the 
two glass slides. For the pure liquid crystals no sol- 
vent was used. For DSC measurements, the samples 
were prepared by introducing approximately 20-40 
mg of the solution into aluminum DSC pans and the 
samples were left for 24 hours to evaporate the sol- 
vent. Then the pans were covered with an aluminum 
top. 

2.3 DSC Measurements 

DSC measurements were performed on a Perkin 
Elmer Pyris 1 apparatus equipped with a Perkin 
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Elmer 2P intracooler. A heating ramp of 5°C/min 
was used in the temperature range spanning from 10 
to 100°C. Two runs were performed per sample, and 
the data was analyzed in the second run. 



2.4 Optical 
merits 



Microscope Measure- 



A Lcica DMLB microscope equipped with a video 
camera (Leica DC 100) and a hot stage (Linkam 
THMS 600) was used. Samples were heated above 
the solubilization temperature, leading to a single 
isotropic phase, and left for approximately 20 min- 
utes in the homogeneous isotropic state. Afterwards, 
the measurements were performed as follows: 

• For the isotropic-isotropic equilibrium temper- 
ature determination, the microscope was used 
without polarizers (TOM); the sample was 
cooled at three different cooling rates (between 
0.1 and 5°C/min, depending on the sample) until 
the appearance of a second phase was observed 
at each cooling rate. The equilibrium tempera- 
ture was determined by extrapolating to cooling 
rate zero. 



For the isotropic/smectic-A transition tempera- 
ture, cross-polarizers were used (POM), follow- 
ing the same general procedure as before heating 
from the liquid crystalline state until the disap- 
pearance of the smectic phase was observed. 



• In the quenching experiments, the samples 
were cooled at the fastest rate possible (about 
20°C/min), using liquid nitrogen as a coolant. 
Two types of quenches were performed: one-step 
(direct quench from homogeneous to smectic- 
isotropic temperature region) , and two step (one 
quench to the isotropic-isotropic region, hold for 
10 minutes, and then quench to isotropic-smectic 
region) . 



3 Modeling and Simulation 

3.1 Thermodynamic Model 

The theoretical approach used to computationally de- 
termine the smectic PDLC phase diagrams follows 
that presented by Benmouna et al [17] utilizing the 
Flory-Huggins theory of isotropic mixtures [24] ex- 
tended to account for the presence of smectic liq- 
uid crystal phase-ordering using the Maier-Saupe- 
McMillan theory [25]. 



3.1.1 Flory-Huggins theory 

The Flory-Huggins isotropic mixing free energy per 
mole of cell is [24]: 
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where R is the universal gas constant, r^c/^p is the 
ratio of the volume of liquid crystal/polymer with re- 
spect to the reference volume (taken as the volume of 
the smaller species present, in this case a repetitive 
unit of the polymer), \ IS th° mixing interaction pa- 
rameter, (f> is the volumetric fraction of liquid crystal 
(LC), and T is the temperature. The interaction pa- 
rameter is considered a function of temperature only, 
given by \ = A + B/T, and can be determined by 
fitting the experimentally determined Tj/j values. 

3.1.2 Maier-Saupe-McMillan theory 

The equilibrium free energy contribution from smec- 
tic ordering according to Maier-Saupe-McMillan the- 
ory is [25] : 
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where v = A.5AT NI /T (see eqn. 17 of ref. 25), s/a 
is the nematic/smectic order parameter, and £ is the 
smectic interaction parameter. The value of this pa- 
rameter determines the ratio Tnj/Tai, and the order 
(first or second) of nematic/smectic transition [25]. Z 
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is the partition function and is given by [25]: 
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model has already lead to successful predictions of ex- 
perimentally observed phenomena in polymer smec- 
togens [29] and is expected to contribute new insights 
into the structure of smectic PDLCs. 



(6) 

The total free energy density of the mixture is the 
sum of eqns. [4] and [5] 

3.1.3 Phase diagram computation 

For a bi-phasic equilibrium at constant temperature, 
the equilibrium condition is given by equality of the 
chemical potentials of each component in both phases 
(represented by a and (3). 
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In the smectic phase the order parameters must 
minimize the free energy, thus assuming a is the 
smectic phase: 
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The triple point is found as the intersection be- 
tween the I/I and I/SmA curves 

3.2 Material Transformation Model 

The modeling and simulation approach used to com- 
putationally determine nano-scale structure of ex- 
perimentally observed smectic morphologies uses a 
phcnomenological Landau-de Gennes model [27] first 
proposed for the direct isotropic/smectic-A transition 
by Mukherjee, Pleiner, and Brand [28]. This novel 



This model is applicable to pure component 
smectic-A liquid crystals and is used in this work 
based upon a set of experimentally supported as- 
sumptions (given in Section 4.3). This spatial model 
accounts for smectic-A ordering and elasticity and 
does not take into account polymer/liquid crystal in- 
teractions, as in the thermodynamic model (Section 
3.1). An extension to this material transformation 



Equilibrium exists between either two isotropic 
phases (I/I), one isotropic and one smectic (I/SmA), 
or two isotropic and one smectic (triple point). In 
the isotropic phases, a = s = 0. The convergence of 
numerical methods when solving eq. [7] usually fails 
when approaching the critical point. This point can 
be determined using the following condition: 



model to account polymer/liquid crystal interactions 
has not been formulated and is beyond the scope of 
this paper, but as shown later, the pure smectogenic 
model reveals the origin of structural features driven 
by smectic layering. 

3.2.1 Order Parameters 

Theoretical characterization of orientational and 
translational order of the smectic-A mesophase (see 
[TJ) is accomplished using order parameters that ade- 
quately capture the physics involved. Partial orien- 
tational order of the nematic phase is characterized 
using a symmetric traceless quadrupolar tensor [27]: 
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where n/m/1 are the eigenvectors of Q, which char- 
acterize the average molecular orientational axes, and 
S/P are scalars which characterize the extent to 
which the molecules conform to the average orien- 
tational axes [30-32]. 

The one-dimensional translational order of the 
smectic-A mesophase in addition to the orientational 
order found in nematics is characterized through the 
use of primary (orientational) and secondary (trans- 
lational) order parameters together [33]. A complex 
order parameter can be used to characterize transla- 
tional order [27] : 

ipe^ = A + iB (11) 



where <j> is the phase, ip is the scalar amplitude of the 
density modulation, and A/B is the real/imaginary 
component of the complex order parameter. 
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3.2.2 High-Order Landau/de Gennes Model 

A dual-order parameter Landau-de Gennes model for 
the first order isotropic/smectic- A phase transition 
is used that was initially presented by Mukherjee, 
Pleiner, and Brand [27, 28] and later extended by 
adding nematic elastic terms [34,35]: 



where /j, n / ' /i s is the rotational/smectic viscosity, and 
V the control volume. A higher order functional 
derivative must be used due to the second-derivative 
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where / is the free energy density, /o is the free en- 
ergy density of the isotropic phase, terms 1-5 are 
the bulk contributions to the free energy, terms 6- 
7 are couplings of nematic and smectic order; both 
the bulk order and coupling of the nematic director 
and smectic density- wave vector, respectively. Terms 
8-10/11-12 are the nematic/smectic elastic contribu- 
tions to the free energy. T is temperature, T/vj/Ta/ 
are the hypothetical second order transition temper- 
atures for isotropic/nematic and isotropic/smectic- A 
mesophase transitions (refer to ref. 36 for more de- 
tail), and the remaining constants are phenomenolog- 
ical parameters. This free energy density expression 
is a real- valued function of the complex order param- 
eter iff and its complex conjugate V?*, which makes it 
convenient to reformulate using the real and imagi- 
nary parts of the complex order parameter (see eqn 

n). 

The previously mentioned dynamic timescales are 
included in this model via the Landau-Ginzburg 
time-dependent formulation [37]. The general form 
of the time-dependent formulation is as follows [37] : 



3.2.3 Simulation Conditions 

Two-dimensional simulation of the model was per- 
formed using the Galerkin finite element method 
(Comsol Multiphysics). Quadratic Lagrange basis 
(^Junctions W ere used for the Q-tensor components and 
quartic Hermite basis functions used for the com- 
plex order parameter components. Standard numer- 
ical techniques were utilized to ensure convergence 
and stability of the solution including an adaptive 
backward-difference formula implicit time integration 
method. A triangular mesh was used such that there 
was a density of 6 nodes per equilibrium smectic-A 
layer in the circular computational domain ([2]). 

In order to simulate a boundary condition corre- 
sponding to a bulk texture, surface/boundary inter- 
action with the liquid crystal was modeled assuming 
liquid crystalline ordering at the boundary fully de- 
termined by bulk conditions: 
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where f s (J/m 2 ) is the surface free energy density 
which is coupled to the bulk using using an elastic 
assumption, where the timescale of surface ordering 
is much greater that bulk ordering and the surface 
interactions are assumed always at equilibrium [37]: 
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where k is the unit normal to the surface (outward) 
and 9 is the order parameter. Due to the use of 
the full complex order parameter in the high-order 
model eqn |12| boundary conditions in rp cannot be 
directly used. Instead, boundary conditions must be 
formulated in terms of A/B, the real and imaginary 
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Boundary 



Figure 2: Schematic of the two-dimensional simu- 
lation domain where the initial smectic-A nucleus 
growth direction into the unstable isotropic phase 
is shown. Dark solid lines indicate surface bound- 
aries, grey-shaded regions indicate initially smectic- 
A ordered area, and white indicates initially disor- 
dered/isotropic area. The model parameters used 
in the simulations are as follows: Tmi — 322.85-ff, 
T AI = 330.5X, a Q = 2 x 10 5 ^, b = 2.823 x 10 7 ^ , 
c= 1.972 xl0 7 -^,a 



1.903 xlO 6 - 



3.956 x 



..»K ' @ ' 

e = 1.938 x I0~ n P N, b x = 



10 8 ^3 A = 9.792 xlO 6 ^ 
\x\Q- 12 pN,b 2 = 3.334x10~ 30 Jto, h = lxlO~ 12 pN 
and l 2 = 1.033 x 10~ 12 piV. The ratio of the rotational 
and diffusional viscosities used was — = 25. 



parts of the full complex order parameter. The re- 
sulting bo und ary conditions used in simulation in- 
clude eqn 18 where 9 = A, — 0, B = 0, and 
These boundary conditions correspond 



V 2 S = 

to the assumption that the smectic-A phase <f> at 
the boundary is fixed/constant. Referring back to 
the complex order parameter expression, this corre- 
sponds to specifying the smectic phase <f> = (2n+ 1)tt 
which fixes A = ip and B — at that point. 

In this work full time-dependent simulation of the 
confined smectic domain and initial conditions shown 
in [2] is used to obtain stationary solutions to the 
nano-scalc liquid crystalline structure. The pseudo- 
steady-statc nano-scalc structures are the main focus 
of these simulations and thus the dynamic growth 
and formation mechanisms are not discussed. This 



approach to stationary structure determination has 
a key advantage over relaxation techniques [38] in 
that all stationary solutions are obtained via a phys- 
ically consistent pathway as predicted by the model. 
The main assumptions in this approach are the initial 
conditions (liquid crystal nucleation at the confined 
surfaces) and that thermal effects and flow do not 
play a significant role. 



4 Results and Discussion 

4.1 Phase diagrams 

The experimental (symbols) and theoretical (lines) 
phase diagrams for PS-10CB and PS-12CB, mea- 
sured by POM, TOM and DSC are shown in [3}i,b. 
The diagram for both mesogens shows a homoge- 
neous isotropic phase (I) , coexistence of two isotropic 
phases (I/I), coexistence of a smectic and an isotropic 
phase (I/SmA), and coexistence of a crystalline and 
an isotropic phase (I/Cr). Isotropic/isotropic (I/I) 
phase separation was determined using TOM, the 
isotropic/smectic-A (I/SmA) transition by POM and 
DSC, and crystallization by DSC only. Adequate 
agreement was found between the measurements by 
POM and DSC for the isotropic/smectic-A transi- 
tion, but the DSC results showed a large dispersion, 
so only POM results are shown (in [3]) . For liquid 
crystal concentrations of 50% in weight and greater, 
isotropic/isotropic phase equilibrium was observed. 
In this range of concentrations, the isotropic/smectic- 
A transition temperature and the crystallization tem- 
perature of the liquid crystal-rich phase are com- 
parable to the pure liquid crystal values, indicating 
that the phase separated LC-rich domains are essen- 
tially pure liquid crystal. The increased molecular 
weight of the liquid crystal molecule (12CB versus 
10CB) shifts the phase diagram to higher tempera- 
tures where both Tai and T/j increase. 

|3]also shows the computationally determined phase 
equilibrium results (full and dashed lines) using the 
model and experimentally determined phase diagram 
data. The value of T N/I for 10CB and 12CB (note 
that it is a theoretical second order isotropic/nematic 
temperature [36]) was obtained from literature val- 
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Figure 3: Experimental (points) and theoretical (lines) phase diagrams of (a) PS-10CB and (b) PS-12CB. The 
different regions identified in the diagrams are: homogeneous isotropic phase (I), coexistence of two isotropic 
phases (I/I), coexistence of a smectic and an isotropic phase (I/SmA), and coexistence of a crystalline and 
an isotropic phase (I/Cr). The transition from an homogeneous isotropic state to the I/I region (■) was 
measured by TOM; the transition from I/I to I/SmA region (a) and from an homogeneous isotropic state 
to the I/SmA region (►) were measured by POM and DSC, and the melting temperature (♦) was measured 
by DSC. 



ues [36]. Smectic interaction parameters £ for both 
10CB and 12CB were determined to agree with the 
pure component isotropic/smectic-A transition tem- 
perature ([3j. The mixing interaction parameters 
X = A + B /T were determined via least squares fit- 
ting of the experimentally determined Tj/j and are 
shown in [T] The errors in the interaction param- 
eters were obtained according to error propagation 
theory [39,40]. The computed phase diagram results 
show good agreement with the experimental data. 

Experimental and calculated phase diagrams are 
similar to those reported by Benmouna et al. [16, 
17] for PS/8CB (octyl-cyanobiphenyl) , but with two 
important differences: 

1. From the experimental point of view, 
their system showed the sequence 
isotropic/nematic/smectic, while the present 
system shows a direct isotropic/smectic transi- 
tion. 

2. From the modeling point of view, they have three 
unknown, adjustable parameters A, B and r p , 



while in this work an a priori value for r p is used, 
calculated from its definition, consequently we 
have only two adjustable parameters. 

The fitting method presently used also different. Ben- 
mouna et al [17] use the experimental critical point 
data to fit two of the parameters where in this work 
the two parameters are fit to all data points simul- 
taneously. An interesting finding, as can be seen in 
[l] is that the interaction parameter can be consid- 
ered independent of the molecular weight of the LC 
(within statistical uncertainty) , which is typically not 
the case for low molar mass species. 

Based upon the approach of Benmouna et al [16], 
the DSC data (shown in |4| was utilized to estimate 
the solubility limit of the liquid crystal in the polymer 
[3 (defined as the maximum amount of liquid crys- 
tal that can be dissolved in the polymer at Taj, ex- 
pressed as a weight fraction) and the relative amount 
of liquid crystal in the liquid crystal-rich domains a 
(equal to the ratio of the weights of liquid crystal 
in the liquid crystal-rich domains over the total liq- 
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Table 1: Material properties and model parameters 



parameter 


10CB 


12CB 


M ref 


104-^ 

7TLOL 


104 

711, Oi 


Pref 


1.047^3 

cm" 


1.047-^3 

cm" 


M-e f 


99.33^ 

mol 


99.33^ 

mol 


M PS 


20800 -S 

mol 


20800^ 


PPS 


1.047^3 

cm A 


1-047^ 


Vps 


19860^ 

mol 


19860^ 


M LC 


319^4 

mol 


347^4 
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uid crystal). The solubility limit can be determined 
from the change of enthalpy at the isotropic/smectic- 
A transition for the smectic PDLC material and that 
of the pure liquid crystal [16]: 



AH?P LC 



(19) 



where 5 is the ratio of the enthalpies, H^P LC is the 
enthalpy of phase-ordering of the PDLC material, 
and H%j is the enthalpy of phase-ordering of the 
pure liquid crystal. This ratio 5 can be related to 
the solubility limit (3 using the rule of inverse seg- 
ments [16]: 

where wlc is the weight fraction of liquid crystal in 
the PDLC material. This approach is valid given the 
assumptions [16] that the thermodynamic properties 
of the liquid crystal in both domains (liquid crystal- 
and polymer-rich) are the same, the amount of liquid 
crystal dissolved in the polymer is constant for all ma- 
terials where wlc > 0i & n d the isotropic/smectic- A 
transition temperature is not dependent on wlc (well 
supported by the observations presented in pk,b). 



[5j;,d show 5 versus wlc for both smectic PDLC ma- 
terials (determined using eqn. |19| , where extrapo- 
lation to 8 = yields the solubility limit fi (at a 
reference temperature of Tai from [3^,,b) . The values 
of = 0.40 for PS/12CB and = 0.42 for PS/10CB 
were determined. The fraction of liquid crystal in 
both polymer- and liquid crystal-rich domains can 
also be calculated from 5 [16], under the assumption 
that the liquid crystal-rich phase is essentially pure: 



(21) 



'>LC 



where r is the mass fraction liquid crystal in the LC- 
rich phase over the total liquid crystal mass in the 
material. This can also be calculated again using the 
rule of inverse segments [16]: 



He -13 
He (i - 0) 



(22) 



f5fe,f show r versus He for both smectic PDLC mate- 



rial (determined using eqn. 21 1 qualitatively agreeing 
well with past results [16] where r approaches zero 
as the solubility limit is approached (and no LC-rich 
phase exists) and reaches one as the polymer content 
approaches zero. 

4.2 Macro-scale Morphologies 

The determined phase diagrams ([3| provide quali- 
tative guidance on which phase will be continuous 
and which will be dispersed, or if both will be co- 
continuous (depending if the concentration is lower, 
higher or very close to the critical concentration) , but 
is not sufficient to make quantitative predictions on 
the actual morphologies. In the case of PDLCs, the 
presence of both phase separation and phase ordering 
transitions results in several different dynamic path- 
ways to affect domain morphology of the final mate- 
rial. Three main categories of factors are involved: 
phase transition sequence, phase transition kinetics, 
and the introduction of anisotropy/texture from the 
presence of the liquid crystal. Phase transition se- 
quence, in the context of the experimentally observed 
equilibrium phase behavior of smectic PDLCs (|3| , in- 
volves either a staggered initiation of phase separa- 
tion/phase ordering (when Tjj and Tai exist) or a 
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(a) (b) (c) (d) 

Figure 6: POM images of LC/PDLC morphologies obtained by cooling from 52°C to 50°C at 0.5°C/min: 
(a) 100.0%wi 10CB, (b) 95.59W 10CB, (c) 75%wt 10CB, (d) 49.79W 10CB. 
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Figure 4: Plots of DSC data for (a) PS/10CB mix- 
tures (from bottom to top 100.0%, 75.0%, 49.7%, 
30.0%) and (b) PS/12CB mixtures (from bottom to 
top 100%, 75.6%, 49.5%, 31.1%). 



simultaneous initiation (when Tj/ai exists, or when 
the sample is rapidly cooled such that I/I phase sep- 
aration is negligible). The presence of anisotropy in 
these materials, as mentioned before, introduces an 
additional degree of complexity. 

A unique attribute to liquid crystals exhibiting 
a direct disordered/ordered smectic transition are 
the complex set of growth morphologies observed, 
which have been the focus of much past study [32] 
(Figure [HJ;). The smectic- A mesophase, with its 
lamellar ordering on the molecular scale, exhibits 
growth, defect, and texture phenomena not seen 
in the nematic phases, including the lamellar-like 
cholesteric mesophase. In addition to spheroids 
and sphero-cylinders observed in ncmatics, direct 
isotropic/smectic-A liquid crystals exhibit the fas- 
cinating highly-anisotropic "batonnet" structures of 
Friedel and Grandjean [41-44] and complex buckling 
filaments [45-49]. As there are many factors involved, 
a discussion of the phase transformation processes in 
these materials will be presented before analyzing the 
experimental results in detail, in order to guide the 
interpretation of the observed morphologies. 

For a smectic PDLC material where both Tjj and 
Tai exist, quenching from the fully isotropic one- 
phase state below Tjj results in phase separation 
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Figure 5: Plots of the enthalpy ratios 
percent liquid crystal for (a) PS/10CB, 
and plots of mass fraction of LC in 
(smectic-A) domains versus total mass 
for (c) PS/10CB, (d) PS/12CB. Curve 



correspond to eqn 20 and for (c)-(d) 
eqn. fusing f3 = 0.42 for PS/10CB 



versus weight 
(b) PS/12CB, 
phase-ordered 
fraction of LC 
fits in (a)-(b) 
correspond to 
and (3 = 0.40 



for PS/12CB). 



occurring into fully isotropic liquid crystal-rich and 
polymer-rich domains. These domains formed during 
phase separation into two different isotropic phases 
result in what will be referred to as "stage I" PDLC 
morphology. In this stage spherical domains are 
formed, unaffected by whether the polymer-rich or 
liquid crystal-rich phase is the matrix phase, in that 
both phases are isotropic. This process is well un- 
derstood as a standard liquid/liquid demixing pro- 
cess where the size and number density of the spher- 
ical domains formed are determined by liquid/liquid 
interactions and quench rate [50-52]. As domains 
formed in stage I are quenched below Tai, their mor- 
phology is typically unaffected by the phase-ordering 
transition and standard smectic transition kinetics 
and coarsening are observed [53] (within the liquid 
crystal-rich domain). In the case where the liquid 
crystal-rich phase is the matrix phase (concentration 
of LC higher than the critical concentration, see^), 
the spherical polymer-rich domains serve to prohibit 
the formation of a neat liquid crystalline texture and 
the final texture can be predicted given that the liq- 
uid crystal/polymer interface preferred anchoring is 
known [6,54,55]. 

Two different phase transition sequences can re- 
sult in stage I spherical morphologies: a) quench- 
ing a fully isotropic PDLC below Tjj held for long 
times (complete phase separation) and b) quench- 
ing a fully isotropic PDLC below Tjj held for short 
times and then quenching again below Tai- When 
the material is cooled below Tai, the liquid crystal 
diffuses from already formed isotropic polymer-rich 
domains to neighboring liquid crystal-rich domains, 
but if this diffusion is slow compared to the cooling 
rate, the isotropic polymer-rich phase becomes super- 
saturated in the LC component. This results in sec- 
ondary phase separation, with formation of new do- 
mains, which are referred to as "stage II," and include 
both newly formed domains (post-stage I) and small 
domains formed during stage I (lengths on the order 
of the characteristic length of liquid crystalline order- 
ing) . Phase separation in this regime results in a myr- 
iad of complex morphologies in that surface tension 
anisotropy is present, in addition to the surface elas- 
tic effects of bulk liquid crystalline texture. Much ex- 
perimental research has been done on characterizing 
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growth morphologies of both pure and lyotropic liq- 
uid crystals which exhibit a direct isotropic/smectic- 
A transition [19, 45-47, 56, 57]. Much of the in- 
teresting morphologies exhibited, including spheri- 
cal, sphero-cylindrical, and filament domains [45-47] , 
are attributed to the presence of interfacial tension 
anisotropy A7 = 7|| — 7j_ which is coupled to the 
concentration of solvent present (or polymer in the 
case of PDLCs) and degree of liquid crystallinity. 

The final type of morphological growth occurs at 
longer time-scales than stage I-II growth, where liq- 
uid crystal-rich spherulites in the micron/sub- micron 
range form in polymer-rich domains that are ap- 
proaching equilibrium concentration. In this case, 
diffusion of the liquid crystal is increasingly less en- 
ergetic, becoming the kinetically controlling step and 
limiting the liquid crystal-rich domain size that can 
be locally formed. As diffusion becomes predomi- 
nant in this regime, the domains formed are not only 
smaller, but also more symmetric than the stage II 
domains. This stage is favored at high polymer con- 
centration, where the molecular mobility is lower. 

For PDLC materials with polymer concentrations 
higher than the maximum solubility limit (as de- 
fined previously) both phase separation and order- 
ing occur simultaneous and Tj/aj is found (refer to 
|3j. In this case, stage I morphologies are not ob- 
served in that isotropic/isotropic phase separation 
does not take place. Thus depending on the relative 
ratio of polymer to liquid crystal, stage II morpholo- 
gies are observed (favored in mixtures with relatively 
less polymer) in addition to stage III morphologies 
(favored in mixtures with relatively more polymer). 
PDLCs in these concentration ranges are more dif- 
ficult to study experimentally because of the small 
domain sizes formed (order of (j,m). 

A simple scaling argument can be derived to de- 
scribe the characteristic micron size of a drop formed 
in stage III. The simplest LC shape equation is a bal- 
ance of the bulk stress jump and anisotropic surface 
tension [32]: 

<9k 2 



(k ■ T b ) ■ k 



7ls 



ei 



(23) 



where — (k • Tb) • k is the stress jump, k is the unit 
normal, 7 is the interfacial tension, H the average cur- 



vature, D the deviatoric curvature , I s = eiei +e2e2 
is the unit surface tensor, q = eiei — e2e2 is the 
deviatoric unit tensor, and e^ are the principal eigen- 
vectors of the curvature tensor ei . For spherical 
shapes ( D = 0, H = — l/R) and neglecting anchor- 
ing effects (^J : I = 0) the drop radius R is 

27 



R = 



(k ■ T b ) ■ k 



(24) 



The interfacial tension 7 and bulk stress jump 
(k • Tb) • k can be expressed in terms of the homoge- 
neous free energy density fh'- 



7 : 



2y/L^ f [A(*,c)-/h(* e , 



c e )] 1/2 (25) 



where the subscript e denotes the bulk equilibrium 
value, and is the gradient elasticity for the smectic 
order parameter. Replacing these expressions into 
equation [24] we find that the spherulite radius R is 
given by a balance between the bulk energy with the 
gradient energy: 



R — C \ 



kT\c* 



C 



J° e 2[/,QM-A(* e ,Ce)] 1/2 <M 
/h(*e,C e ) 



(26) 

where the scaled free energy is fh = fh/kT\c* — c e \ 
and c* is the concentration at which the isotropic 
phase is unstable. Using C = 1, we find that at 
coexistence R is given by : 



(27) 




Assuming w lpN and kT\c* — c e \ w 1-^ eqn. 27 
provides an estimate of R sa \pan. Once the droplet 
forms, the growth process is controlled by diffusion. 
Under diffusion- limited conditions, further growth is 
inhibited and the spherulites remain in the low /im 
range (see[6}i and [7 



4.2.1 Quasi-equilibrium morphologies: 
centration effects 



con- 



[6]ind [7] show representative quasi-equilibrium mor- 
phologies obtained for mixtures with different con- 
centrations of 10CB and 12CB, respectively. The 
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(d) 



Figure 7: POM images of LC/PDLC morphologies obtained by cooling from 60°C to 55°C at 0.5°C/min. 
(a) 1009W 12CB, (b) 95.3%wt 12CB, (c) 62.5%wt 12CB, (d) 30%wt 12CB. 



samples were heated from room temperature to a 
temperature above the isotropic/smectic-A transition 
(60°C for mixtures with 12CB and 52°C for mixtures 
with 10CB), and then they were cooled at 0.5°C/min 
and held at the final temperature for 20 minutes be- 
fore the images were taken. The initial temperatures 
were below the isotropic/isotropic solubilization tem- 
perature for mixtures with more than 50%wt liquid 
crystal, such that some degree of phase separation 
took place before quench below the phase-ordering 
transition. [6^ and [7^, show the morphologies of the 
pure liquid crystal, where the typical focal conic tex- 
ture of a smectic phase can be observed [41]. As poly- 
mer concentration is increased, a two-phase morphol- 
ogy can be seen. [6]d and|7jD show the morphologies of 
samples with 95%wt LC (above the the critical-point 
concentration, which was calculated to be 0.87 for 
10CB and 0.88 for 12CB) where the liquid crystal- 
rich phase is the matrix. Isotropic/isotropic phase 
separation produces a droplet morphology (stage I) 
where the spherical dispersed polymer-rich domains 
form in the liquid crystal-rich matrix. Once quenched 
below Taj, the continuous liquid crystal-rich matrix 
phase transitions to smectic-A via the standard nu- 
cleation/growth mechanism with additional topolog- 
ical constraints imposed by the dispersed polymer- 



rich domains (see previous discussion). Addition- 
ally, small liquid crystal-rich domains are observed 
to form inside the isotropic polymer-rich droplets, 
giving rise to a "salami"-type structure frequently 
observed in polymer mixtures (for example, high- 
impact polystyrene [58,59]). The formation of these 
morphologies is due to secondary phase separation, 
as explained before (stages II-III). 



When the concentration of liquid crystal becomes 
lower than the critical concentration, the isotropic 
phase is the matrix and the smectic phase is dis- 
persed. The size of the smectic domains becomes 
smaller as the concentration of liquid crystal de- 
creases. In this sub-critical concentration range of 
liquid crystal, the morphology of the material is com- 
posed of a combination of large and small domains, 
where the large domains are predominantly spher- 
ical in shape and the small domains have a much 
more asymmetric shape. The spherical domains were 
formed by the I/I phase separation (stage I), and the 
smaller and asymmetric domains correspond to the 
stages II-III. In[6jl and[7]i, the mixture composition 
is such that there is no I/I equilibrium, so stage I 
morphologies are not observed. 
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4.2.2 Effects of quench rate on morphology 
dynamics 

In order to further study the formation and dynamics 
of different morphologies, two types of rapid quenches 
from a temperature above 7/j were performed: (a) 
direct quench at 20°C/min to the isotropic/smectic- 
A coexistence region and (b) quench at 20°C/min to 
the isotropic/isotropic region, hold for 10 minutes, 
quench at 20°C/min to the isotropic/smectic-A re- 
gion. In quench (a) isotropic/isotropic phase sepa- 
ration is minimized such that stage I morphologies 
are rare, whereas in quench (b) isotropic/isotropic 
phase-separation is allowed to progress, resulting in 
many spherical stage I domains. [8] shows the time- 
progression of the domain morphologies at different 
times during the quench sequence (a) (directly into 
the isotropic/smectic-A regime), for a mixture con- 
taining 85%wt of 12CB. It is observed that, in this 
case, the vast majority of smectic-A domains formed 
correspond to stage II morphology. [8] shows the time- 
progression of quench (b) where distinct stages I- 
III are formed, directly exhibiting the influence of 
isotropic/isotropic phase separation preceding phase- 
ordering. In the direct quench, when the phase sepa- 
ration/transition process begins the matrix has a rel- 
atively high concentration of liquid crystal, such that 
a large majority of stage II domains form. In the 
two-stage quench, due to isotropic/isotropic phase- 
separation, secondary phase separation/transition 
takes place after the final quench in a matrix phase 
with a relatively high concentration of polymer, such 
that stage III morphologies are prevalent also. 

4.3 Nano-scale Structure 

In the final stage of the study, the material trans- 
formation model was simulated in a two-dimensional 
circular smectic-A domain. Initial and boundary 



conditions (see Sections 3.2 1 were used to simulate 



a bulk radial-oriented smectic-A domain, similar to 
past studies [38,48]. Previously determined phe- 
nomenological parameters were used based on 12CB 
[60] to perform a parametric study in temperature, 
with fixed cylindrical length scale of 97.5nm. Time- 
dependent simulation was performed until quasi- 



equilibrium textures were formed, after simulation 
times of > I [is. [9] shows the three different sta- 
tionary smectic-A textures predicted by the model. 
Two general nano-scale structures are predicted, ei- 
ther a liquid crystalline defect core ([9^, and ^p) or 
an isotropic/disordered core p]:) depending on the 
degree of undercooling. 

The fully smectic cores ([9^ and [9J5) are similar to 
the bipolar radial textures seen in confined nemat- 
ics [31], where a pair of +| disclinations form at the 
core to mediate the LC topology imposed by the com- 
bination of intcrfacial confinement geometry and an- 
choring. These orientational defects have biaxial core 
structures [61], which can be identified through the 
degree of biaxiality [62,63]: 



B 2 



1 - 6 



Kg ■ Q) ■ QY 
(Q ■ Q? 



(28) 



which ranges from — ► 1 (no biaxiality to maxi- 
mum biaxiality). [9^i shows the smectic layer struc- 
ture that is predicted by the model for intermedi- 
ate and deep undercooling. The core is composed 
of an eccentric cylindrical/continuous smectic layer 
with the pair of + | disclination (orientational) de- 
fects along the major axis. These texture and de- 
fect structures are consistent with past work study- 
ing the transient growth of an initially radial-textured 
smectic-A spherulite [60]. Past results of cylindrical 
confinement using a lower order model [38] predicted 
the presence of a smectic +1 disclination, which was 
not observed to be stable under the conditions in this 
work. 

|9]d shows the smectic layer structure as tempera- 
ture is increased, where expansion of the smectic layer 
spacing collapses the cylindrical core smectic layer 
into a misoriented layer that undergoes a high degree 
of layer expansion at its center. [T0^ an d [Top show 
the scalar smectic order parameter eqn. 1 1 1 1 for both 
fully smectic core structures. These results agree 
with past work in nematic systems where an increase 
in temperature increases the separation distance be- 
tween the disclination defects [31]. In addition to de- 
fects in the underlying orientational order, the model 
predicts the formation of translational defects, ele- 
mentary dislocations, enclosing the misoriented core 
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(d) (e) (f) 



Figure 8: POM images of LC/PDLC morphologies obtained after quenching a mixture with 85.1% 12CB. 
starting at 80°C, quench at 20°C/min to 50°C. the snapshots were taken at (a) 1:25, (b) 2:45, and (c) 10:00 
minutes. Quench starting at 80°C, quench at 20°C/min to 53°C, hold for 10 minutes, quench at 20°C/min 
to 50°C. The snapshots were taken at (d) 0:14, (e) 1:46, and (f) 10:00 minutes 




(a) (b) (c) 



Figure 9: Simulation results of the smectic density modulation Re("i>) (surface) and biaxiality (contour) 1 : 
(a) T = 331.00JT, (b) T = 331.35JC, and (c) T = 331.41X. Surface plot values range from Re{9) = -0.161 
(black) to Re(^) = 0.161 (white) and the contour corresponds to j3 2 = 0.1. 
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layer. These observations elucidate structural differ- 
ences, on the nano-scale, of cylindrical smectic do- 
mains. 

The second general cylindrical domain structure 
([9^, and 10 1) shows a fully melted/disordered core 
structure that has been predicted in the past for 
smectic-A filaments [48]. In this shallow undercool- 
ing scenario, the free energy difference between the 
smectic and isotropic phases is less than that of the 
formation of defects, as in the cases of [9^-b. 

The two-dimensional simulation results pre- 
sented also provide approximations of fully three- 
dimensional spheroidal and sphero-cylindrical do- 
main structure. These three-dimensional geometries 
share symmetries with the presently simulated two- 
dimensional geometry and thus through projections 
on and rotations around these symmetry axes full 
three-dimensional nano-scale structure can be in- 
ferred. QT] shows schematics of rotational axes of 
symmetry between a two-dimensional circle (present 
simulation geometry) and the extrusion/rotation re- 
lations with cylinders, spheres, and sphero-cylindcrs. 

The first relation of the two-dimensional domain 
[TTk and [TTb to a three-dimensional filament [TTb in- 
volves a projection into the z-axis. This relation 
takes advantage of the two-dimensional simulation 
assumption of gradients in the z-axis equaling zero. 
Using the nano-scale structure predicted in [9^,, a di- 
rect quantitative filament structure is predicted to be 
composed of a eccentric cylindrical smectic-A layer 
capped by +1/2 disclination lines. The structure 
predicted at intermediate temperatures [9|d results in 
a disoriented smectic layer capped by both disclina- 
tions and elementary dislocation lines. 

The second relation of the two-dimensional do- 
main |11^ to a three-dimensional spheroid |llb in- 
volves a rotation on the x/y-axis. This relation 
provides only qualitative information on the pos- 
sible three-dimensional spherical domain structure. 
A direct qualitative relationship would require two- 
dimensional axial-symmetric simulation where the 
model is reformulated in spherical coordinates and 
subsequent gradients in the azimuthal axis are as- 
sumed zero. Past work simulating nematic mate- 
rials using both two-dimensional rectangular (the 
present approach) and axial-symmetric formulations 



has shown qualitative structural agreement [64,65], 
supporting the present assumption. Thus, these ro- 
tations result in two possible smectic layer structures 
with either point or line defects (high-energy point 
defect structures are neglected [65]). Using the nano- 
scale structure predicted in [9^, the resulting spheroid 
structure is composed of a prolate spherical smectic- 
A layer with a single circular +1/2 disclination loop 
about its major axis. The structure predicted at in- 
termediate temperatures results in a disoriented 
disc-shaped smectic-A layer confined by single circu- 
lar disclination and elementary dislocation loops. 

Finally the sphero-cylindrical structure 11: which 
results from a combination of the structures inferred 
from relations shown in |llb and |llp l involves capping 
smectic-A filaments with half-spherical domains. In 
these cases, defect loops are elliptical instead of cir- 
cular. 



5 Conclusions 

In conclusion, a combined experimen- 
tal/computational study of phase equilibrium 
and growth morphologies of novel PDLC mate- 
rials which exhibit a direct isotropic/smectic-A 
transition was performed. Partial phase diagrams 
were determined for two different smectic PDLC 
materials: PS/10CB and PS/12CB. The approach of 
Benmouna et al [17] utilizing the Flory-Huggins the- 
ory of isotropic mixing and Maier-Saupe-McMillan 
theory for smectic-A liquid crystalline ordering was 
used to computationally determine phase diagrams 
for both systems. The resulting computationally 
determined phase diagrams were found to be in good 
agreement with the experimental results, and provide 
both an estimation of the full phase diagram and the 
ability to predict the effects of using mono-disperse 
PS of different molar masses [17]. 

In addition to thermodynamic data, specific trends 
of the phase transition morphologies and dynam- 
ics were found for different compositions/phase- 
transition sequences and a general theory in- 
troduced identifying three distinct morphological 
regimes. These regimes include the formation 
of essentially spherical domains (isotropic/isotropic 
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(a) (b) (c) 



(d) 

Figure 10: Simulation results of the smectic scalar order parameter tp (surface) and biaxiality (contour) 1 : 
(a) T = 331. 00K, (b) T = 331.35K, and (c) T = 331A1K. Surface plot values range from tp = -> 0.161 
corresponding to (d) the color bar and the contour corresponds to 1 = 0.1. 



phase separation, "stage 1"), highly anisotropic do- 
mains (isotropic/smectic-A phase separation, "stage 
2" ) , and sub- micron size spheroid domains (late 
isotropic/smectic-A phase separation, "stage 3"). 

Finally, two-dimensional simulation of a high-order 
Landau-de Gennes model was used to determine 
nano-scalc structure of generalized stage 1-2 morpolo- 
gies. A parametric study in temperature was per- 
formed under three different degrees of undercool- 
ing. These novel results for two-dimensional smec- 
tic spherulites are consistent with past work in nc- 
matics [31], where the core structure is composed 
of +1/2 disclination loops. In addition to agreeing 
with past results, the presented simulations predict 
the smectic-A layer structure of the core to both in- 
clude (intermediate undercooling) and lack (deep un- 
dercooling) elementary dislocation loops. Symmetry 
relations were then used to relate two-dimensional 
simulation results to fully three-dimensional exper- 
imentally observed morphologies, predicting possi- 
ble nano-scale structures of spheroidal and sphero- 
cylindrical domains. 

This study provides a basis for future experimental 
and theoretical study of these novel smectic PDLC 



materials which exhibit a direct isotropic/smectic- 
A phase transition. The next step of experimental 
analysis (using the techniques of TOM/POM) in- 
volve increased resolution studies of quench dynamics 
and growth morphologies, guided by the presented 
phase diagrams (|3^,b). In addition to experimen- 
tal approaches, there is great promise in modeling 
and simulation of these systems, where length and 
time scales (nano-scale) inaccessible experimentally 
are able to be resolved. In particular, the extension 
of presently used Landau-de Gennes type model of 
Mukherjee, Pleiner, and Brand [28] to take into ac- 
count the presence of polymers would be a substantial 
first step. 
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(d) (e) (f) 



Figure 11: Schematics of (a) two-dimensional simulation results for intermediate/deep quenches [9p, (x/y/z- 
axes correspond to horizontal/ vertical/out-of-plane axes) (b)-(c) the z-axis projection of the two-dimensional 
structure to a full three-dimensional cylindrical/filament domain (d)-(e) the y-axis rotation of the two- 
dimensional structure to a full three-dimensional spherical domain (f) sphero-cylindrical combination of 
z-axis projection (c) and y-axis rotation (e). 
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